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Abstract. Let X = {X(t),t > 0} be a Brownian motion or a spectrally negative 
stable process of index 1 < a < 2. Let E = {E(t),t > 0} be the hitting time of 
a stable subordinator of index < p < 1 independent of X. We use a connection 
between X(E(t)) and the stable subordinator of index P/a to derive information 
on the path behavior of X(E t ). This is an extension of the connection of iterated 
Brownian motion and (i)-stable subordinator due to Bertoin [JJ. Using this con- 
nection, we obtain various laws of the iterated logarithm for X(E(t)). In particular, 
we establish law of the iterated logarithm for local time Brownian motion, X(L(t)), 
where A is a Brownian motion (the case a = 2) and L(t) is the local time at zero of 
a stable process Y of index 1 < 7 < 2 independent of X. In this case E(pt) = L(t) 
with P = 1 — I/7 for some constant p > 0. This establishes the lower bound in the 
law of the iterated logarithm which we could not prove with the techniques of our 
paper [57] ■ We also obtain exact small ball probability for X(E t ) using ideas from 



Key words and phrases. Local time Brownian motion, hitting time of stable subordinator, spec- 
trally negative stable process, law of the iterated logarithm, Hirsch's integral test, Kolmogorov's 
integral test, small ball probability. 
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1. Introduction 



Self-similar processes arise naturally in limit theorems of random walks and other 
stochastic processes, and they have been applied to model various phenomena in a 
wide range of scientific areas including telecommunications, turbulence, image pro- 
cessing and finance. In this paper, we will study a class of self similar processes that 
are non-Gaussian and non-Markovian. 

To define the processes studied in this paper, let X = {X(t),t > 0} be a Brownian 
motion or a spectrally negative stable Levy process of index 1 < a < 2 started 
at 0. When a = 2, we denote X = W. Let E = {E(t),t > 0} be the hitting 
time of a stable subordinator of index < (3 < 1 independent of X. Our aim in 
this paper is to obtain various laws of the iterated logarithm (LIL) for the process 
Z = {Z(t) = X(E(t)), t > 0}. This process is self-similar with index H = f3/a. 

Many authors have constructed and investigated various classes of non-Gaussian 
self-similar processes. See, for example, [38] for information on self-similar stable 
processes with stationary increments. Let W be an independent copy of W. Define 
B\t) = W(t),t > and B\t) = W'(-t), t < 0. Burdzy [131 E] introduced 
the so-called iterated Brownian motion (IBM), A(t) = B 1 (B(t)), by replacing the 
time parameter in B l (t) by an independent one- dimensional Brownian motion B = 
{B(t),t > 0}. His work inspired many researchers to explore the connections between 
IBM (or other iterated processes) and PDEs, to establish potential theoretical results 
and to study its sample path properties, see [H [3H [32J, [3U [351 EH1 HQ] and references 
therein. 

The process Z emerges as the scaling limit of a continuous time random walk 
with heavy-tailed waiting times between jumps [51 [29]. Moreover, Z has a close 
connection to fractional partial differential equations. Baeumer and Meerschaert [3], 
and Meerschaert and Sheffler [2S] showed that the process Z can be applied to pro- 
vide a solution to the fractional Cauchy problem. More precisely, they proved that, 
u(t,x) = E, x [f(Z(t))} solves the following fractional in time PDE 

(1.1) —u(t,x) = G x u(t,x); u(0,x) = f(x), 

where G x is the generator of the semigroup of X t and d^gfy/dt 13 is the Caputo 
fractional derivative in time, which can be defined as the inverse Laplace transform 
of s^gi^s) — s /3 ~ 1 (7(0), where g(s) = J °° e~ st g{t)dt is the usual Laplace transform. Let 
L(t) be the local time at zero of a strictly stable process Y(t) of index 1 < 7 < 2 
independent of X. The process, W(L(t)) is called local time Brownian motion in [27]. 
Recently Baeumer, Meerschaert and Nane [1] further established the equivalence of 
the governing PDEs of X(L(t)) and X(\B(t)\) when (3 = 1/2. Here B = {B(t),t > 0} 
is another Brownian motion independent of X and Y. When 7 = 2, the process 
W(L(t)) has also appeared in the works of Borodin [101 [11], Ikeda and Watanabe 
[23] . Kasahara [21], and Papanicolaou et al. [37] . 
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We are interested in sample path properties of Z. Since Z is non-Gaussian, non- 
Markovian and does not have stationary increments, the existing theories on Markov 
and/or Gaussian processes can not be applied to Z directly and some new tools need 
to be developed. In this regard, we mention that Csaki, Foldes and Revesz [18] 
studied the Strassen type LIL of Z(t) = W(L(t)) when L(t) is the local time at zero 
of a symmetric stable Levy process; Meerschaert, Nane and Xiao [27] established a 
large deviation result, and uniform modulus of continuity for W(L(t)). In particular 
they obtained upper bounds in the LIL and the modulus of continuity of W(L(t)). 
In this paper, we will establish LIL results for Z. 

This paper is a continuation of [TS] and [27J. In section 2, we study small ball 
probability for Z(t). Our results in Section 2 imply in particular that there can not 
be a Chung type LIL for Z due to power decay in the small ball probability. In 
section 3, we study asymptotic properties of Z — {Z(t), t > 0} by first establishing 
a connection to a stable subordinator of index (3/oc, which is an extension of the 
connection of iterated Brownian motion and ^-stable subordinator due to Bertoin [7J 
and then making use of limiting theorems on stable subordinators. 



The study of small ball probability is useful for studying Chung type LIL. In this 
section we show that there cannot be a Chung type LIL for Z. 

In what follows we will write f(t) ~ g(t) as t — > a to mean that for some positive 
G\ and C2, C\ < liminfi^ a f(t)/g(t) < limsup t _, f(t)/g(t) < G%. We will also write 
fit) ~ g(t), as t -> a, to mean lim t _ a f(t)/g(t) -> 1. 

Next we introduce the processes that will be studied in this paper. A Levy process 
X = {X(t),t > 0} with values in R is called strictly stable of index a G (0, 2], a ^ 1 
if its characteristic function is given by 



where — 1 < v < 1 and x > are constants. In the terminology of [2SJ Definition 
1.1.6], —v and are respectively the skewness and scale parameters of the stable 

random variable X(l). When a = 2 and x = 2, X is a standard Brownian motion 
and, if x = 1, X is a Brownian motion running at twice the speed of standard 
Brownian motion and we denote it by W. When 1 < a < 2 and v — 1, X is called 
spectrally negative stable Levy process. 

In general, many properties of stable Levy processes can be characterized by the 
parameters a, v and \. For a systematic account on Levy processes we refer to [6]. 

Let E(t) be the continuous inverse of a stable subordinator D t of index < (3 < 1: 



2. Small ball Probability 



(2.1) 




E(t) 



mi{s : D(s) > t}, (t > 0). 



3 



E(t) has also been called the hitting time of stable subordinator of index < (3 < 1 
in the literature, see, for example, [5]. 

The next theorem is our first result on the small ball probability of W(E(t)). 

Theorem 2.1. Let W(t) be a Brownian motion. Let E(t) be the continuous inverse 
of a stable subordinator D{t) of index < (3 < 1. Then 



limu^P sup \W(E(t))\ < u 



0<t<l 



32r( / 3) ski^vr) 



71 ' 



\k-i 



j(2*-l) 



Proof. From a well-known formula (see Chung pL6] ) : 



4 

sup W(t) < w ] = - 

0<*<1 7T 



k=l 



(2k - 1) 2 7T 2 



Since -E(t) is continuous and nondecreasing, we use conditioning to get 



sup \Z(t)\ < u 

0<i<l 



sup \W(t)\<u\ 

0<KE(1) 



E 



P sup \W(t)\< 

o<t<i y E(l) 



E{1) 



(2.2) 



fc-i 



7T — 2k- 1 

fc=i 



■ exp 



(2fc - l) 2 7r 2 £(l 



The Laplace transform of (-E'(i) inverse stable subordinator of index (3 ) (cf. 
Remark 2. 8]) is 



(2.3) 



E[eMsE(t))) = = J2 f^^) 



n=0 



which is the Mittag-Leffler function evaluated at — st 13 . And, by equation (13)], 
for some c > 



(2.4) 



< E p {-at p ) <c/(l + ab- 



using Equation (12.41) and the fact that Mittag-Leffler function is completely monotone 
we obtain that the series in (I2.2p is absolutely convergent, we get by dominated 
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convergence 



E 



~ /_ 1)fc -i 

k=l 



E 



fc-i 



E( exp 



7T ^ 2A; - 1 V 

fc=i v 



\fc-i 



7T ^ 2Jfe - 1 
fc=l 




Therefore, using the fact that 



(2.5) 
with c 



z ~ cz as z 



oo 



r(/3) sin(/37r) 



see, Krageloh [26] and taking the limit w 



using dominated convergence we obtain 
suPo< f <i \Z{t)\ < u ] 



lim ■ 

u->0 



4 ~ (_i)fc-i 

— > urn 

7T ^— ' 2A; - 1 n^O 

fc=l 



■ inside the sum 

(2fc-l) 2 7T 2 ^ 



8u 2 



:-ir 



fc-i 



(2.6) 



vr^ 2k -1 (2k-l) 2 n 2 
32c 

(2fc- l) 3 



□ 



Remark 2.2. Since the exact small ball probability in Theorem 12. II is a power decay, 
we do not expect to get a Chung type LIL for Z(t). We will obtain a Hirsch's integral 
test below for Z(t) using another method. 

Using the same conditioning method in Theorem 12.11 and ideas in Aurzada and 
Lifshits [2] and Nane [33], we also can get the following 

Theorem 2.3. Let U = {U(t),t > 0} be a self similar process of index < H < 1, 
and let E{t) be the inverse of a stable subordinator of index < (3 < 1, independent 
ofU. Let 9 = 1/H > and suppose 



(2.7) 
Then 



logP ( sup \U(t)\ < u ) ~ ku~\ as u 
y o<t<i 



0. 



P ( sup \U(E(t))\ < u 

K 0<t<l 



m 1 ^ as w 



0. 



Remark 2.4. In this theorem we get upper and lower bounds for the small ball prob- 
ability due to the lack of an explicit formula for the small ball probability for the 
process U. 

Proof. Let u > 0. By assumption, for all < e < e = €o(v), 

-k(l+v)u- e <■ p / \TT(f\\ <■ < -k(i-v)u-° 



< P sup \U(t)\ < u < e 
\o<t<i / 

This implies that there are constants ci,C2 > that depend on eo such that for all 
u > 0, 

(2.8) Cl e~ k{l+U)U ' B < P ( sup \U(t)\ < u j < C2e-*( 1_,/)u ~ fl 

\o<t<i / 

Since is continuous and nondecreasing, we use conditioning to get 
P[ sup \U(E(t))\ < u ] = P[ sup |E/i|<«,] 

0<t<l 0<t<B(l) 



E 



P( sup \U(t)\ < TvTwW) 

0<i<l -^UJ 



(2.9) < c 2 S [exp (-fc(l - z/)M" e E(l))] . 

we use the upper bound in (12. 8p to get the last inequality. We also obtain a lower 
bound for the small ball probability of U(E) using the lower bound in (12.81) . 

Now taking limit as u — >• in (12. 9p using equations (12.31) . (12. 5p and then letting 
v — > we get the upper bound for the small ball probability of U(E). With a similar 
argument we obtain the lower bound. □ 

Fractional Brownian motion (fBm) is a centered Gaussian process W H = {W (t), t G 
M} with ^^(O) = and covariance function 



(2.10) E(W H (s)W H (t)) = ^{\s\ 2H + \t\ 



2H \„ +\2H 



where H e (0, 1) is a constant. W H is a self-similar process of index H. The small 
ball probability is given for W H with 9 = 1/H in (JM]), see, for example, [38J. We 
obtain the small ball probability for local time fractional Brownian motion, W H (L(t)), 
here L(t) is the Local time at zero of an independent strictly stable process of index 
1 < 7 < 2, as a consequence of the following corollary using the fact that E(pt) = L(t) 
with j3 = 1 — I/7 for some constant p > 0, see, for example, [39] . 

Corollary 2.5. Let W H = {W H (t),t > 0} be a fractional Brownian motion with 
index < H < 1, and let E(t) be the inverse of a stable subordinator of index 
< P < 1, independent ofW H . Then 

P ( sup \W H {E(t))\ < u j w as u — ► 0. 
\o<t<i / 
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The process W H (E) emerges as the scaling limit of a correlated continuous time 
random walk with heavy-tailed waiting times between jumps |28j . 

Strictly stable process of index < a < 2 is a self-similar process of index \ja 
and the small ball probability is given with 9 = a in (12.71) . see, for example, [7J. We 
obtain the small ball probability estimate of Z as a consequence of 

Corollary 2.6. Let X be a strictly stable process of index < a < 2 and let E t be 

the inverse of a stable subordinator of index < (3 < 1, independent of X . Then 

P ( sup \X(E(t))\ <u) ~ u a as u — > 0. 

\0<i<l / 

Remark 2.7. Since the Laplace transform of the process E(t) decays like E(e~ zE ^) ~ 
c(f3)/z as z —>■ oo for which the power of z is independent of the index /?, the small 
ball probabilities does not depend on (3 in the theorems in this section. 

For other extensions of small ball probabilities for iterated (fractional) Brownian 
motions and iterated stable processes, see Aurzada and Lifshits j2]. 

3. Path regularity of Z 

We will consider the process Z = {Z(t), t > 0} defined by 

Z(t) = X(E(t)), Vt>0, 

where X = W for a = 2 and X is the spectrally negative stable Levy process with 
1 < a < 2, v — 1. 

Let Sit) = sup 0<s<t X (s) . Since E{t) is an increasing and continuous process, we 
have 

Z{t) = sup Z(s) = S(E(t)). 

0<s<t 

Let a a be a stable subordinator of index I /a, then 5* is the continuous inverse a a , 
see Bingham [H Proposition 1];. 

S(t) = inf{s : a a (s) > t}, (t > 0) 

more precisely, 

E(exp(-sa a (t))) = exp(-t(s/c!) 1/Q ) (s > 0,t > 0), 

where c\ = y -1 sec(7r — |7ra). 

The inverse D(t) of E(t) is a stable subordinator of index < (3 < 1. Here D has 
Laplace transform, ~E(e~ sD ^) = e~ tsP . According to Bochner [9], Do a a is a stable 
subordinator with index ^, and more precisely, 

/•oo 

log (E[e- sDoCTQ W]) = mt / (e-* - l)x-( 1+/3 / Q ^x = -c7 1/a t// Q 
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by integration by parts, we can show that m = c l 1 ^ a (j3/a)[T(l — /3/a)] 1 . Define 

H = (r(l - B/a)m) a/ ^\a - 3)/ 8. 

Let A = (8 / a) / {1 — 8 / a) = [3 /{a — (3). Plainly, Doa a is the right-continuous inverse 
oiSoE. 

Proposition 3.1. The right- continuous inverse of the supremum of Z(t) is 

m£{t : Z(t) > ■} = o- z (-) } 

where o z is a stable subordinator with Laplace exponent s — > c x s^^ a , 

We can now apply certain results in the literature on stable subordinators to derive 
asymptotic results about Z. 

Our first result on LIL is the following 

Theorem 3.2. Let X be a spectrally negative stable Levy process with 1 < a < 2, 
v = 1, and let E(t) be the inverse to a stable subordinator of index < 8 < 1 
independent of X . Then 

UmSUP ^(logllgV"^/- = W ' ^ a ' S ' 
both as t — > + and t — >• oo. Here c(/3, a) = fj}' x 

When a = 2, the upper bound in this large time LIL result was proved by Meer- 
schaert, et al. [27] for the case E pt = L t is the local time at zero of a strictly stable 
process of index 1 < 7 < 2 and in this case (3 — 1 — I/7. Here p > given by 

p -(i-i/ 7 ) = n -^Y(l + l/7)r(l - lh)x~ l Re{[l + wtan(7T7/2)]- 1/7 }, 

see Stone [39]. Hence we obtain the following LIL for local time Brownian motion 
W(L(t)). 

Corollary 3.3. Let W be a Brownian motion, and let L(t) be the local time at zero 
of a stable process of index 1 < 7 < 2, independent ofW. Then 

limsup -. 1V9 ^ <kL P\ u -TW5 = {pc{(3 = 1 - I/7, 2))^^ a.s. 

&ot/i as t — > + anc? t — >• 00. 

Remark 3.4. Corollary 13.31 is an extension of the results in J27J [IS] . Also the case 
7 = 2 gives the LIL for a version of iterated Brownian motion, which was 

obtained in [T31 [T7J, here B is a Browninan motion independent of W . 

Proof of Theorem \3.B. Applying the LIL for stable subordinators proved by Fristedt 
[T9] (see also Breiman [12]) to D o a a , we have 

liminf — j-— — , - - — = u 1//A = c(8, a) a.s. 
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both as t —>■ + and t —>■ oo. 

Since S a o E(D o cr") = t, we deduce that the probability that 

t < S a o E(cr //3 (log | logt|) {/3 -" )//3 ) 

infinitely often as t — > + (or as t ^ oo) equals 1 for c > c(/3, a) and for 
c < c(P,a). Recall that t — »■ t^ a (log | logt|)^ a- ^' a is an asymptotic inverse of 
t ^r/ /3 (log|logt|)( /3 - a )/' 3 , so 

limsup , jr^—rz — sr = (c(/3, at))^ a.s. 
^/ a (log|logt|)( a -^/ Q 

both as £ — ► + and t — ► oo. 

One can replace S o E(t) with □ 

We next obtain the Kolmogorov's integral test for Z. 

Theorem 3.5. Let X be a spectrally negative stable Levy process with 1 < a < 2, 
v — 1, and let E(t) be the inverse to a stable subordinator of index < (3 < 1 
independent of X . Then 

F[Z(t) > t p/a f{t) infinitely often as t -»• oo] = or 1, 

according as 

It 

f(t)V^UM-m)] a/ia - 0) )j 

converges or diverges. A similar result holds for small times by replacing the integral 
over (1, oo) with the integral over (0, 1). 

Proof. Breiman [12] obtained a refinement of Fristedt's LIL for stable subordinators. 
It states that if if : (0, oo) — ► (0, oo) is a decreasing function, then 

(3.1) F[D o a(t) < t Q// V(t) infinitely often as t -> oo] = or 1 

according as the integral 

converges or diverges. Here A = (3 / '(a — (3). 

Let / : (0, oo) — > (0, oo) be an increasing function, then by first replacing t by 
t& a f(t) in (E3]), we get that 

D o a a {t p/a f{t)) < tf{t) a/ ^{t Pla f{t)) infinitely often as t -> oo 

and then taking S o E of both sides we get the probability that 

t p/a f(t) < S o E(tf(t) a ^<p(t^ a f(t)))t -> oo 

equals or 1 according as I(cp) < oo or = oo. If we now choose if such that 
(f(t^ a f(t)) = f(t)- Q//3 , we then see that 

¥[SoE(t)> t p/a f(t) infinitely often as t -»• oo] = or 1, 



according as I((p) < oo or I{tp) = oo. An easy calculation shows that the latter 
alternative reduces to decide whether the integral (replace (p(t) with fit)" / 13 ) 

converges or diverges. We can replace SoEbyXoE = Z, which gives the Kol- 
mogorov's integral test for the upper functions of Z for large times. A similar result 
holds for small times by replacing the integral over (1, oo) with the integral over (0, 1). 

□ 

We next obtain Hirsch's integral test for Z(t) 

Theorem 3.6. Let X be a spectrally negative stable Levy process 1 < a < 2, u = 1, 
and let E(t) be the inverse to a stable subordinator of index < /3 < 1 independent 
of X . Then 

liminf Z(t)/tp(t) = or oo 

t— »oo 

according as the integral (p(t)t~( 1+ P' a > dt diverges or converges. A similar result 
holds for small times by replacing the integral over (1, oo) with the integral over (0, 1). 

Proof. Recall that if / : (0, oo) — ► (0, oo) is an increasing function, then 

limsup D o a a (t)/f (t) = or oo a.s. 

according as f(t)~P/ a dt converges or diverges. See Theorem 11.2 on page 361 in 
Fristedt [20]. Consider first the case Limsup is equal to 0. Let e > 0. Then, 

D o a a (t) < ef(t) for all t large a.s. 

Now taking S o E of both sides we get 

t<So E(ef(t)) for all t large a.s. 

Using scaling of S o E we get 

t < (ef{t)f la t-P /a S o E(t) for all t large a.s. 

hence 

t 1+ P /a (ef(t))-P /a <SoE(t) for all t large a.s. 

We use a similar argument for the case Limsup is infinite. Hence, setting (p(t) = 
t i+/9/«/(t)-/J/« we g e t 

lim inf S o E(t)/(p(t) = or oo 

according as the integral f?° (p(t)t~^ 1+l3 ^ a,> dt diverges or converges. Again, there is a 
similar result for small times. □ 

We next obtain a modulus of continuity result for the supremum process Z. 
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Theorem 3.7. Let X be a spectrally negative stable Levy process with 1 < a < 2, 
v — 1, and let E(t) be the inverse to a stable subordinator of index < j3 < 1 
independent of X . Then 

lim sup - , , T- —. -T7TT? — 5w — d(a, 8) = (ac 2 (5/a)/fl)' 1 ~' 3 ^ a.s. 
^ere c 2 = c 2 (3/a) = (1 - d/a)(d/a)^ a -^ . 

Proof. Since <r z is a stable subordinator of index /3/a, Theorem 1 in Hawkes [2T] gives 
(here we change the Levy measure of a a so that the Laplace exponent is s — >■ s 13 ^) 

lim inf ^(t + h)-a z (t) _ {a „ m 

where c 2 = ciifijot) = (1 — (3 /a)((3/a) t3 ^ a ~ l3 \ Taking inverses as in Hawkes [21] we 
obtain 

lim sup ug f a n + n]h\v}-l)/ a = d ( a ^) = (ac 2 ((3/a)/(3) {1 ~ f3/a) a.s. 

e-»0+ <t<l,0<h<e ftW a [l0g(l//l)p P)l a 

□ 

Remark 3.8. We can also deduce the analog of the Chung- Erdos-Sirao integral test 
for Z from Theorem 7.2.1 in Mijnheer [30 J. We refer to Meerschaert, Nane and Xiao 
[2T] for the modulus of continuity of W(L(t)), where L(t) is the local time at zero 
of an independent strictly stable process Y of index 1 < 7 < 2. We also refer to 
Khoshnevisan and Lewis [25] for the modulus of continuity of IBM A(t). 

Remark 3.9. Bertoin [7], remarked that similar results holds for n-iterated Brownian 
motions. In this regard we get similar LIL and integral tests for n-iterated E(t) 
processes. We do not know whether this is the case for n-iterated spectrally negative 
Levy processes. 

Acknowledgment Author would like to thank Yimin Xiao for discussion of the 
results in this paper. 
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